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Abstract 

We have derived the effect of steady heat fiux on the rate of chemical reaction based on the line- 
of-centers model using the explicit velocity distribution function of the steady-state Boltzmann 
equation for hard-sphere molecules to second order. It is found that the second-order velocity 
distribution function plays an essential role for the calculation of it. We have also compared our 
result with those from the steady-state Bhatnagar-Gross-Krook(BGK) equation and information 
theory, and found no qualitative differences among them. 
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I. INTRODUCTION 



Chemical reactions in gases has been studied with the aid of gas kinetic theory 0, 0,0,0]. 
Present proposed the line-of-centers model for chemical reaction in gases, which has been 
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accepted standard mode 
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to describe the chemical reaction in gases. 
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100,00,0,0, 
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Under nonequilibrium situations such as gases under the heat conduction and the shear 



flow, their nonequi 
among researchers. 



ibrium effects on the rate of chemical reaction have attracted attention 
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23| However, to derive the effect of steady 



heat flux on the rate of chemical reaction in the line-of-centers model, we need the ex- 
plicit velocity distribution function of the steady-state Boltzmann equation for hard-sphere 
molecules to second order in density and temperature gradient. To our knowledge, nobody 
had derived the explicit velocity distribution function of the Boltzmann equation for hard- 
sphere molecules to second order. Although Burnett determined the second-order pressure 

n 

tensor for the Boltzmann equation y|, he had not derived the explicit second-order velocity 
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This is a result of mathematical 



131 ] Therefore, none has succeeded 



distribution function of the Boltzmann equation, 
difficulties, as was indicated by Fort and Cukrowski. 
to obtain the correct reaction rate of Present's model except for Fort and Cukrowski who 
adopted information theory |26j as the nonequilibrium velocity distribution function to sec- 
ond order. Q] 

Recently, we have derived the explicit velocity distribution function of the steady-state 



251 ] In this letter, we apply 



Boltzmann equation for hard-sphere molecules to second order, 
the explicit velocity distribution function to the calculation of the rate of chemical reaction 
in the line-of-centers model under steady heat conduction. 



II. SIMPLE MODEL FOR CHEMICAL REACTION 



In the early stage of a chemical reaction between monatomic molecules: 



A + A — > products, 



(1) 



the rate of chemical reaction is not affected 



viewpoint of kinetic collision theory [ 



wag 



xy the existence of products. 27( From the 
the rate of chemical reaction can be 



2 



described as 



R 



dw I dvi dQ, ffiga(g), 



(2) 



where v and Vi are the velocities of the molecules, g = |v — Vi| their relative speed, O the 
solid angle, / = /(r, v) and f\ = /(r, v x ) are the distributions of v and v x at r, respectively. 



We adopt the line-of-centers mode 

^■flflflflQflyyyy^ 



as the differential cross-section of chemical reaction 
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<?(g) 
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g < 



2l\ The model contains 



I \ mg z I & — y m 



(3) 



with m mass of the molecules and E* the threshold energy of the chemical reaction, d is 
regarded as a distaoee between centers of — .Cedes at contact.QQ 

In order to calculate the rate of chemical reaction (j2J), we expand the velocity distribution 
function / to second order as 



/ = / (o) + / (i) + / ( 2 ) = / (o )(1 + (i) +0 ( 2))) 



(4) 



around the local Maxwellian, = n(m/2iiK,T) 3 / 2 exp[— mv 2 /2k,T], with n the density of 
molecules, k the Boltzmann constant and T the temperature defined from the kinetic energy. 
Substitution of eq.(JH) into eq.(J2J) leads to 



(5) 



up to second order. The zeroth-order term of R 
R(f>) - I ,h, I ,hr. i ,/o / f(o)f(°) 



rfv y rf Vl y <*n y /^i u Vto) = 4nV (™ T 
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e «r. 



(6) 



corresponds to the rate of chemical reaction of the equilibrium theory. Similarly, the first- 
order term of R is obtained as 



RM = JdvJdv.JdnJ /^V + tfVfo). 



(7) 



where R^ does not appear because <ffi^ is an odd functions of c, as will be shown in the 
next section. The second-order term of R, i.e. R^ 2 \ is divided into 



R {2A) = fdvfdvijdnf /^VMVG?), 



(8) 
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and 

R(2 ' B) = I dV I dVl I dn I /( ° )/l(0) ^ (2) + *?Vfo)- (9) 
Since the integrations (JSJ) and © have the cutoff from eq.fjHJ), the explicit forms of 
and of the steady-state Boltzmann equation for hard-sphere molecules are required to 
calculate RP 1 '^ and RS 2,B \ respectively. 



III. EXPLICIT FORM OF THE VELOCITY DISTRIBUTION FUNCTION TO 
SECOND ORDER 

We introduce the explicit form of the velocity distribution function of the steady-state 
Boltzmann equation for hard-sphere molecules to second order in the density and the tem- 



perature gradient along x-axis 



25j: 



^-^tri^f^^r^), (10) 



which is an odd function of c, and 

4096m J 2 „ „, 3, 



(2) 



5Q25b 2 n n 2 K 3 T 3 



{J2r%rT(r + -)Sl(c 2 ) + 



r>2 



E r! M2^. -4- 4)T(r + V f (c 2 )}, (11) 

r>0 

with the scaled velocity c = (m/2KT) 1 / 2 v. Here T(X) represents the Gamma function and 
Sl(X) is a Sonine polynomial, the definition and the orthogonality properties of which are 



written in our paper. 25] J x = — 75&iiK 3 / 2 T 1 / 2 <9 E T/647r 1//2 m 1//2 <i 2 denotes the steady heat flux 
with d x = d/dx. The numerical values for b± r , &o and &2r are listed in Table HI Note that, 
though we have derived them in forms of fractions, we use them in forms of four significant 
figures, since the forms of the fractions are too complicated. We have derived each of <jP^ and 

to 7th approximation of Sonine polynomials, though the complete forms of <fft> and 
for hard-sphere molecules are the sums of an infinite series of Sonine polynomials. However, 
we have confirmed that (p^ does not converge to 4th Sonine approximation but both of 



and (j)^ almost converge to 7th Sonine approximation j25j, so that, in this letter, we will 
show the results calculated from cj)^ and for 7th approximation of Sonine polynomials. 
We note that Burnett had determined only b^r and b^ r to 4th Sonine approximation: 6 2 r is 



related to the second-order pressure tensor. 



24j Therefore, we could not calculate the effect 



of steady heat flux on the rate of chemical reaction with the line-of-centers model if we did 
not derive the explicit forms of and (p^ 2 \ especially &o in ref. 25]. 

In order to compare the results from the steady-state Boltzmann equation with those 
from the steady-state Bhatnagar-Gross-Krook(BGK) equation and information theory, we 
also use the explicit forms of and < - 2 - ) for them. The expressions of and (p^ for 



the steady-state BGK equation can be reduced from the general form o 



28, 



the Chapman- 
while those 



Enskog solution for the steady-state BGK equation to arbitrary order 
for information theory were used by Fort and Cukrowski. Q| We mention that (jP^ for 
the steady-state BGK equation is an odd function of c and identical with that for information 
theory, while are different from each other. 



IV. RESULTS: RATE OF CHEMICAL REACTION 

Inserting and (j)^ 2 ' of eqs.(|10jl and and the corresponding velocity distribution 
functions for the steady-state BGK equation and information theory into eqs.(JBJ) and Q, and 
performing the integrations with the chemical reaction cross-section (j3J) , we finally obtain the 
nonequilibrium effects on the rate of chemical reaction based on the line-of-centers model. 
The expressions of RP ,A ^ and R^ 2 ' B ^ become 

(§)'>■ 

- r>0 V ' 



and 



«,^(^)vs {I >(§y } , 

respectively. The numerical values for a r and (3 r are listed in Tables ITT1 and HTT1 respectively. 
We have found that R( 2 ' B ^ for the steady-state Boltzmann equation is determined only by 
the terms of b 0r in dx_J of eq. (HH) which Burnett [24] had not derived: the terms of &o we 
have derived in ref. [25j are indispensable for the calculation of R^ 2 ' B ^ . 

The explicit expressions in eqs.(|12j) and (|13|) for information theory were already given 
by Fort and CukrowskL_ 12 j_ 13 1 Though they were interested in forms of Qa = R^ 2,A ^ / 
and Q B = R^/R^^M, 

in the present letter, we focus on the nonequilibrium effects 
on the rate of chemical reaction in the forms of R^ 2,A ^ and R^ 2 ' B ^ . This is because the forms 
of Qa and Q b for the steady-state Boltzmann equation do not converge even when we adopt 
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any higher order approximation of Sonine polynomials. The upper limit of r in eqs. ()12|) and 
(fTHj) is directly related to the order of the approximation of Sonine polynomials in eqs. (jlO|) 
and (fTTj): ^(c 2 ) includes the term of c 2r . We emphasize that, however, the values of RS 2 < A "> 
and R( 2,B ^ for the steady-state Boltzmann equation converge to 7th approximation of Sonine 
polynomials. It should be also mentioned that there appear artificial qualitative differences 
in both Qa and Qb for the steady-state Boltzmann equation, the steady-state BGK equation 
and information theory for large E*/kT. 

The graphical results of compared with those of R^ 2 ' A ^ are provided in Fig^ Both 
of R {2) and R (2 > A) in Figfflare scaled by vr 1 / 2 rf 2 m 1 / 2 J 2 /K 5 / 2 T 5 / 2 . Note that R {2) is the sum of 
R( 2 ' A "> and i?( 2 ' B ) in eqs. (fT2|) and (JT3J). As FigHJshows, it is clear that R^ 2,B ^ plays an essential 
role for the evaluation of R^ 2 \ We have also found that there are no qualitative differences in 
R^ of the steady-state Boltzmann equation, the steady-state BGK equation and information 
theory, while R^ 2 ' A ^ also exhibits a slight deviation from each other. This deviation in 
r( 2 > a ) would not be observed if we adopted (ffi-' of the steady-state Boltzmann equation for 
the lowest approximation of Sonine polynomials as in the previous papers jl^l ll^ . This is 
because (f)^ 1 ' of the steady-state Boltzmann equation for the lowest Sonine approximation is 
identical with the precise ft 1 ' of the steady-state BGK equation and information theory. It 
should be mentioned that we have found qualitative differences among these theories in some 



301 ] The nonequilibrium 



thermodynamic quantities, e.g. the second-order pressure tensor. |25l. 
effects on the rate of chemical reaction is an insensitive quantity to the differences among 
the three theories. 



V. DISCUSSION 

0( 2 ) is indispensable for the calculation of the nonequilibrium effects on the rate of chem- 
ical reaction, since RS 1 ^ does not appear and R( 2 ' B ^ is remarkably larger than R 1 - 2 '^ as FigQ 
shows. Thus, the nonequilibrium effect on the rate of chemical reaction will substantiate 
significance of the second-order coefficients in the solution of the steady-state Boltzmann 
equation, although their importance has been demonstrated only for descriptions of shock 
wave profiles and sound propagation phenomena jsil. We should emphasize that 

the calculation of the second-order pressure tensor and the application of it to shock wave 
profiles and sound propagation do not require the explicit form of the second-order veloc- 
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ity distribution function: the second-order coefficients derived by Burnett do not include 
any terms of &o in <ft^ of eq. (jll|) . 24, |2| The nonequilibrium effects we consider appear 
mainly in the ranges of small E* / kT where the chemical reaction often occurs, so that our 
results should be compared with experimental results in the early stage of the chemical 
reaction)^; a large additional effect due to some modification of the velocity distribution 
function because of the chemical reaction itself becomes significant as the chemical reaction 
proceeds Q,Q- 

We also propose a thermometer of a monatomic dilute gas system under steady heat 
flux. We mean that we can measure the temperature T around a heat bath at To in the 
nonequilibrium steady-state system indirectly with the aid of the nonequilibrium effect on 
the rate of chemical reaction as follows. The nonequilibrium effect on the rate of chemical 
reaction in the early stage around the heat bath can be measured. Thus, one can compare the 
experimental result with the theoretical result shown in FigQwith T = T . The difference 
between the former and the latter will indicate that the temperature T around the heat 
bath is not identical with T , that is, T = T + A. Here A should depend upon the heat flux 
in general. Substituting this temperature expression into the explicit expressions of m 
eq.©, we obtain a new correction term 

as the nonequilibrium effect on the rate of chemical reaction besides in FigQ with 
T = To. At last, we can estimate the gap A so as to make the new correction term match 
the experimental result. For example, if A in eq.([_14p is proportional to the heat flux, 



we will confirm the relevancy of the slip effect [34. l35| . If A in eq. (|14)) is identical with 

2mJ 2 /5n 2 K 3 T 2 , i.e. 1/T = 1/T — 2mJ 2 /5n 2 fi; 3 T 4 , the experimental result will agree with 

the theoretical result from the steady-state Boltzmann equation with the nonequilibrium 

f~| 

temperature 9 = T predicted by Jou et al.[2ty. We show the comparison of the theoretical 
results, i.e. R^ for the steady-state Boltzmann equation with T = To and that with 
6 = T , in FigH where i? (2) is scaled by vr 1 / 2 c/ 2 m 1 / 2 J 2 /k 5 / 2 T 5/2 . We note that one can 
obtain the explicit form of R^ with 9 = T expressed as the dashed line in FigO as the 
sum of R^ from eqs.(|12|) and (|13|) for the steady-state Boltzmann equation with T = T 
and R {new) from eq.flHJ) with A = 2mJ 2 /5n 2 /t 3 Tq. We have found that there is a significant 
difference between R^ with T = T and that with 9 = T for small E*/kT . This significant 



difference also appears in for the steady-state BGK equation and information theory. 
Similar difference in for information theory would also follow from the results of Fort 
and Cukrowski ["K^ . although they were not interested in the ranges of small E*/kT. We 
emphasize that, however, our proposal written in this paragraph differs from that by Fort 
and Cukrowski [13( which uses a given value measured directly by a thermometer in order 
to determine T or 9 regardless of the temperature of the heat bath. 
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(2) 

Scaled R 



R 




Boltzmann Eq. 
BGK Eq. 

information theory 



R 

Boltzmann Eq. 

BGK Eq.and information theory 



E* /kT 



FIG. 1: Scaled compared to scaled R( 2 ' A ^ as a function of E* / kT for the line-of-centers model. 




E*/kT c 



FIG. 2: Comparison of scaled from the steady-state Boltzmann equation with T = T$ (the 
solid line) and that with 9 = Tq (the dashed line) as a function of E*/kTq. Note that the former 
is identical with the result shown by the solid line in Fig.l with T = Tq. 
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TABLE I: Numerical constants bi r , &or and &2r in eas. tfTU)) and 



T 


Plr 


OOr 


02r 


n 
u 


n 
U 


1 
1 


1 ^90 v i n — 2 


i 
i 


1 09^ 


n 
u 


1 97R v 1 D~ 1 


9 
z 


4 8Q9 v 1 D~ 2 


ti.OOU X 1U 


fi in y 1 n — 2 


3 


3.715 x 1CT 3 


-5.429 x 1CT 2 


5.521 x 10~ 3 


4 


2.922 x 10 -4 


-4.098 x 10~ 3 


4.214 x 10~ 4 


5 


2.187 x 10~ 5 


-3.184 x 10" 4 


3.106 x 10~ 5 


6 


1.492 x 10~ 6 


-2.087 x 10~ 5 


1.861 x 10~ 6 


7 


8.322 x 10~ 8 







TABLE II: Numerical constants a r in ea. (|12|) . 



r 


Boltzmann Eq. 


BGK Eq. and Information Theory 





-2.292 x 10~ 2 


l 

50 


1 


-1.448 x lO" 1 


3 

25 


2 


3.223 x 10" 1 


6 

25 


3 


-9.834 x 10 -2 


4 
75 


4 


7.919 x 10" 3 




5 


-6.752 x 10" 4 




6 


5.298 x 10~ 5 




7 


-3.584 x 10~ 6 




8 


2.039 x 10~ 7 




9 


-9.599 x 10~ 9 




10 


3.613 x 10" 10 




11 


-1.032 x lO -11 




12 


2.110 x 10~ 13 




13 


-2.887 x 10~ 15 




14 


2.352 x 10~ 17 




15 


-8.578 x 10" 20 
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TABLE III: Numerical constants r in ea. (fT31) . 



/• 


Boltzmann Eq. 


BGK Eq. 


Information Theory 





-1.361 x 10 _1 


n 

75 


19 
150 


1 


-5.094 x 1CT 1 


38 
75 


7 
15 


2 


3.968 x KT 1 


4 
15 


26 
75 


3 


5.805 x 10~ 2 


8 
75 


4 
75 


4 


-2.811 x 10~ 3 






5 


1.039 x 10" 4 






6 


-1.808 x 10~ 6 
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